Non-Abelian Stokes Theorem and Quark- Mo nop ole 

Interaction!] 



Dmitri Diakonov and Victor Petrov 

Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188350, Russia 

Abstract 

We derive a new non-abelian Stokes theorem by rewriting the Wilson loop as a 
gauge-invariant area integral, at the price of integrating over an auxiliary field from 
the coset SU(N)/[(U(1)] N ~ 1 space. We then introduce the relativistic quark-monopole 
interaction as a Wess-Zumino-type action, and extend it to the non-abelian case. We 
show that condensation of monopoles and confinement can be investigated in terms of 
the behaviour of the monopole world lines. One can thus avoid hard problems of how 
to introduce monopole fields and dual Yang-Mills potentials. 



1 Introduction 

On the way to describe confinement as due to the monopole condensation several 
serious technical problems arise. Until now monopole condensation has been fully 
theoretically understood only in theories undergoing spontaneous breaking of colour 
symmetry down to the £7(1) subgroup(s). We mean a) the Georgi-Glashow model in 
2+1 dimensions by Polyakov Q and b) the supersymmetric model in 3+1 dimensions 
by Seiberg and Witten Q. In both cases there is an elementary Higgs field in the 
adjoint representation, whose non-zero vacuum expectation value mercilessly breakes 
the gauge group down to the U(l). The compactness of the remaining U(l) group 
allows Polyakov-'t Hooft monopoles. In 2+1 dimensions they are 'pseudoparticles' and 
cannot condense; in 3+1 dimensions they can. In both cases a linear confining potential 
is obtained for particles which are electrically charged in respect to the unbroken U{\) 
subgroup. 

x The extended version is published in: Nonperturbative approaches to QCD, Proceedings of the Internat. 
workshop at ECT*, Trento, July 10-29, 1995, D.Diakonov (ed.), PNPI (1995) 



In principle, in QCD a similar scenario could take place: the role of the elementary 
Higgs field breaking the gauge symmetry could be played by some composite gluon 
operator belonging to the adjoint representation; its v.e.v. might break the colour 
SU(3) down to U(l)xU (I). However, such a possibility seems to be ruled out by what 
we know experimentally: we would have far more hadrons than in reality. Therefore, 
one should probably assume that the colour group remains unbroken - in contrast to 
the supersymmetric example of Seiberg and Witten 0. 

An introduction of monopoles without elementary or dynamical Higgs mechanism is 
a difficult task by itself. To override the difficulty, 't Hooft || has suggested to visualize 
the U(l) monopoles by partially fixing the gauge, up to the U(l) transformations. This 
procedure is called abelian projection. Monopoles are then objects similar to those 
which one would find in the symmetry-broken case. Of course, such objects are to a 
great extent dependent on the concrete choice of the gauge, or the concrete choice of 
the abelian projection used. Their desirable condensation is even more obscure: even 
if it happens in one gauge, it need not necessarily happen in another Therefore, it 
would be helpful to introduce a gauge-invariant description of monopoles and a gauge- 
invariant formulation of the monopole condensation, if there is any. This paper is a 
step in that direction. 

It should be noted that, if the colour group remains unbroken, the mere notion 
of "monopole condensation" becomes somewhat ambiguous. It seems, at least, that 
it can not be the usual Higgs mechanism (applied to dual Yang-Mills potentials), for 
the same reasons: there would be a proliferation of hadron states. Indeed, two (dual) 
gluons out of eight would be 'electrically neutral' in respect to the two U(l) subgroups 
where the monopoles live, as well as two new kinds of mesons with the colour structure 
t/SA 3 ' 8 V> and five new kinds of baryons which happen to have the colour T3 and Y 
equal to zero. All these states (multiplied by the variety of quantum numbers) are not 
affected by the dual Meissner effect, therefore they are not confined and should be thus 
observable as new types of hadrons pi. It is difficult to imagine that this new realm of 
hadrons has somehow avoided registration, the more so that, as recently observed in 
ref. ||, an additional discrete symmetry makes some of these unusual hadrons stable 

2 There is an endless philological discussion whether the Higgs mechanism actually breakes gauge symme- 
try. Fortunately, there exists a gauge-invariant formulation, for example: Consider a two-point correlation 
function C(x — y) = (cj) a (x) P a b[x , y)(j) b (y)) , where 4> a is a scalar (possibly composite) field and P a b is a P- 
exponent in the adjoint representation. C(x — y) is perfectly gauge invariant. If C(x — y) decays exponentially 
at large separations, there is no symmetry breaking; if it tends to a constant, the situation is referred to as 
Higgs mechanism and spontaneous gauge symmetry breaking; if it decays as a power of the separation, it 
can witness the Berezinsky-Kosterlitz-Thouless phase. 
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under strong interactions. 

If confinement is not due to the dual Higgs effect, then what is due to? There is in- 
creasing evidence from lattice investigations that monopoles extracted by the maximal 
abelian projection are important in getting the area law for the Wilson loop. Therefore, 
a direct theoretical study of the relation between monopoles and confinement is highly 
desirable, without referring to such evasive quantities as the monopole field and the 
dual Yang-Mills potentials. For that reason we prefer the first-quantization formalism 
for monopoles (i.e. path integrals over monopole loops) rather than the field-theoretical 
one. Trying to avoid the notion of the dual fields, we derive the direct (but of course 
non-local) interaction between non-abelian charges and non-abelian monopoles. 

We start by rewriting the Wilson loop as the ordinary exponent of a certain flux 
through the surface spanned on the closed contour of a heavy quark loop. The price 
one has to pay for that is an additional integration over an auxiliary scalar field n from 
the SU{N)/[U(1)]( N -V coset space. (In the N = 2 case n is a unit 3- vector, n 2 = 1). 
We call it the non-abelian Stokes theorem. It is a smart formula, in fact. In order to 
fix the representation to which the probe quark of the Wilson loop belongs, we have to 
add a Wess-Zumino-type term for the n field. With that term, the flux becomes that 
of a gauge-invariant field strength introduced earlier by Polyakov |7| and 't Hooft ||] in 
relation to monopoles, and the Stokes theorem applies to that particular field strength. 
The auxiliary n field plays the role of the direction of the elementary Higgs field in 
colour space. 

We next present a relativistic formula for the interaction of a point charge and 
a point monopole. It is also a Wess-Zumino-type formula, but in three dimensions. 
Using the results of our previous work j|] we are able to formulate the quark-monopole 
interaction for the non-abelian case. 

This paper deals with " kinematical" problems, leaving aside the hard dynamical 
one: what is the actual driving force for confinement. However, we hope that this paper 
gives a framework for a gauge-invariant description of the confinement mechanism 
based on monopole condensation, whatever it means. 
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2 Non-abelian Stokes theorem 



The path ordering in the Wilson loop, 



L = Tr P exp \ i j A^dx^j , (1) 
can be eliminated at the price of introducing integration over all gauge transformations 



along the loop [10|. Let r parametrize the loop and A{t) be the tangent component 
of the Yang-Mills field along the loop in the fundamental representation of the gauge 
group, A(t) = A^dx^/dT, Tr(t a t b ) = \5 ah . The gauge transformation of A(t) is 



d 



A{t) S(t)A(t)S~\t) + iS(r) — S~\t). 

ar 



(2) 



Let Hi be the generators from the Cartan subalgebra (i = 1, ...,r; r is the rank of the 
gauge group) and the r-dimensional vector m be the eldest weight of the representation 
in which the Wilson loop is considered. The formula for the Wilson loop derived in ref. 
[ p~C| ] is a path integral over all gauge transformations S(t) which should be periodic 
along the contour: 



L = J DS{t) exp ijdrTr mi Hi (SAS' 1 + iSS' 
For example, in the simple case of the SU (2) group eq. (||) reads: 



L = / DS(t) exp 



U J dT Tr t 3 {SAS^ + iS&) 



(3) 



(4) 



where r 3 is the third Pauli matrix and J = t> , 1, |, ... is the 'spin' of the representation 
of the Wilson loop considered. In what follows we shall concentrate for simplicity on 
the SU(2) gauge group. 

The path integrals over all gauge rotations (^, ||) are not of the Feynman type: 
they do not contain terms quadratic in the derivatives in r. Therefore, a certain regu- 
larization is understood in these equations: for example, one can introduce quadratic 
terms in the angular velocities iSS^ with small coefficients which eventually should be 
put to zero; see ref. [10 1 for details, where eq. (||) have been derived in two independent 
ways. Let us stress that eqs.(|3],[|) are manifestly gauge invariant, as is the Wilson loop 
itself. 

The second term in the exponent of eqs.(QQ) is in fact a Wess-Zumino-type term, 
and it can be artificially rewritten not as a line but as a surface integral inside the 
closed contour of the Wilson loop. Indeed, let us parametrize the SU{2) matrix S 
from eq. (|j) by the Euler angles, 
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S = exp(iVr 3 /2) exp(i9r 2 /2) exp(^r 3 /2). (5) 
The second term in the exponent of eq. (0) is then 

i J /" dr Tr^SS^) = iJ / dr(cos 9 <j) + ip) (6) 

where the last term is in fact zero due to the periodicity of S(t). 
Introducing a unit three vector 

1 . „ 

n = -Tr (StS't ) = (sin 9 cos 4>, sin 9 sin 0, cos 9) (7) 
we can rewrite (||) as 

(6) =i~J dTdae abc e ij n a d l n b d j n c , (8) 

where one integrates over any surface spanned on the contour. Indeed, the integrand 
of eq. @ is known to be a full derivative; using the Stokes theorem (the standard one!) 
one reproduces eq. ©. Let us note that the r.h.s. of eq. (|8|) is the 'topological charge' 
of the n field at the surface: 

Q = — [ dTdae abc e ij n a d i n b d j n c . (9) 

Eq. (^) can be also rewritten in a form which is invariant under the reparametriza- 
tions of the surface. Introducing the invariant element of a surface, 

e*V = dadr - = ^ dHArea), (10) 

one can rewrite eq. (||) as 

(6) = i~y cPa^n^dpn . (11) 
We get for the Wilson loop 



L = J Dn(a,r) exp iJ J dT{A a n a ) + y / d 2 a a pe abc n a d a n b d l 



in ' 



(12) 



The interpretation of this formula is obvious: the unit vector n plays the role of 
the instant direction of the colour "spin" in colour space; however, multiplying its 
length by J does not yet guarantee that we deal with a true quantum state from a 
representation labelled by J - that is achieved only by introducing the Wess-Zumino 
term in eq. ([l2|): it fixes the representation to which the probe quark of the Wilson 
loop belongs to be exactly J. 
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Finally, we can rewrite the exponent in this formula so that both terms appear to 
be surface integrals: 



L = J Dn(a,r) exp y (- J d 2 a aP F^n a + J d 2 a af3 e abc n a (D a n) b (D^n)^ , (13) 

where D ab = d a 5 ab + e acb A c a is the covariant derivative and = d a A® - dp A a v + 
e abc ^6 j g ^Yxq ggj^ strength. Indeed, expanding the exponent of eq. ( |13| ) in powers 
of A a one observes that the quadratic term cancels out while the linear one is a full 
derivative reproducing the A a n a term in eq. (|l2|); the zero-order term is the Wess- 
Zumino term (||). Note that both terms in eq. ([l3|) are explicitly gauge invariant. We 
call eq. (^) the new non-abelian Stokes theorem. Another version of a non-abelian 



Stokes theorem has been suggested some time ago by Simonov [13]. 
One can introduce a gauge-invariant field strength, 

G a p = F^n a - e abc n a (D a n) b {D p n) c , (14) 

which, as a matter of fact, coincides in form with the gauge-invariant field strength 
introduced by Polyakov Q and 't Hooft ||] in connection with monopoles. In that case 
the unit-vector field n had the meaning of the direction of the elementary Higgs field, 



3 Quark— monopole interaction 

We start with considering a non-relativistic abelian electric charge e moving in the 
field of a magnetic monopole sitting at the origin and having the magnetic field B = 
(g/47r)r/|r| 3 . The equation of motion for the charge is given by the Lorentz force: 

m Xl = —eyfc-j-p (15) 

It is known that the Dirac quantization condition requires eg = Ann; we shall choose 
n = 1. The relativistic generalization of eq. (|l5|), when the charge is moving along 
some world line x^{t\) and the monopole is moving along some world line ^(ts), is 

d 2 f (x(n) - y(T 2 )) a dxp{n) cfy 7 (r 2 ) 

dTiyJxf, ir J |x(t-i) — y(-r 2 )| dri dr 2 

Indeed, taking the monopole sitting at the origin so that its world line is y^fo) = 
(t2, 0,0,0), and a non-relativistic charge with x M (ri) = (ti, Xj(ri)), \xi\ <C 1, we 
return, after integrating over r 2 from minus to plus infinity, to eq. (|1|). We note that 
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eq. ([!(]) is Lorentz invariant and also invariant under the re-parametrization of both 
world lines, x^{j\) and y^T-z), as it should be. 

The l.h.s. of eq. ( |I6| ) is obviously the variation of the standard relativistic action 
of a free particle, 



Sfree = m J dT V%^ ( 17 ) 

what about the right-hand side? The appropriate interaction term appears to be a 
tricky thing. It can be written only as a non-uniquely defined Wess-Zumino action. 
Let us introduce a unit 4-vector 

U fl {T 1 ,T 2 ) = —— —-, U =1, 18 

\%{n) - y{T 2 )\ 

and define its analytical continuation to a third dimension labelled by a, < a < 1, 
so that 



Vp(Ti,T 2 ,(T = 1) = « M (ri,r 2 ), v li (ri,T 2 ,a = 0) = w M , (19) 

where u/,, is some constant 4-vector of unit length. One can also introduce an SU(2) 
unitary matrix 

V(n,T2,a) = Vp(Ti,T2,a) ■ a~, a~ = (1 2 , (20) 

where are the three Pauli matrices and 1 2 is the unity 2x2 matrix, and define three 
anti-hermitean matrices, 

L A = V*d A V, A = Tl ,T 2 ,a. (21) 

The needed relativistic charge-monopole interaction term can be written as the 
winding number of V (times 2n, to make exp(iSi n t) uniquely defined): 

2tt 



24vr 2 



J dTidT 2 daeABC r R (L A L B L C ) 



67T J dT ^ dT ^ dae ABC( aPlS dAVadBVpdcV^VS- (22) 



Indeed, let us vary this S{ n t in respect to the trajectory of the charge x^(ti). To that 
end we first find the change of (p!q ) due to the arbitrary variation 5v^(ti, t 2 , a). We 
have 

eABC^ 5 dAVadsvpdcv^vs = Seabc e a ^^ &a 3v a Ob Vf3 dc v-y vs 

+ eABC^^dAVadBVpdcv-fSvg. (23) 
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The last term is zero here since partial derivatives as well as the variation of u« must 
be all orthogonal to and to each other (because of the antisymmetric e ^ 6 ), which 
is impossible in four dimensions. Therefore, we are left with the first term in eq. (p3|), 
which, for the same reason, can be written as a full derivative, 

(23) = 38 A [e A Bce afhS 6v a d B v l 3dcK r vs] . (24) 

Integrals over full derivatives in ri 2 are zero, if we assume that the trajectories Zu(ri) 
and y M (r2) are closed, so that we are left only with the full derivative in the auxiliary 
dimension labelled by a. Therefore we get 

SSi n t = -~ J dridT2 J da ~^ e^^Vad^vpdr^vs . (25) 

The value of the square brackets at a = is zero since we have chosen the continuation 
( |j~9| ) in such a way that at a = is equal to a constant vector and is thus t\^- 
independent. Therefore, the integral over a reduces to the value of the square bracket 
at a = 1 where assumes its physical value u^(ti, t 2 ), as defined by eq. fllS)). We have 
thus demonstrated a well-known general fact that, though the Wess-Zumino term is 
not uniquely defined, its variation is. 

We now take the variation 5v^ (equal to 5u^ at the physical surface a = 1) to be 
due to the variation of the charge trajectory x^(ti). In the environment of eq. (^) it 
means that 



We get finally 



Rn )-y{T 2 )\ 



SSint 2 f (g(n) - y(T 2 )) a dxpjn) dy 7 (r 2 ) 



Sx^n) ir J \x(n) - y(r 2 )\ dr\ dr 2 

that is the r.h.s of the equation of motion (|l6|), q.e.d. Actually, eq. ( |22"| ) is a general- 
ization of the non-relativistic charge-monopole interaction, see, e.g., refs.[|ll], 12]. 



The main point is that one cannot write down the action whose variation is given 
by eq. (^) , in a unique way. This is because it is impossible to write the electric charge 
- magnetic charge interaction without introducing a string (or other more complicated 
objects) which take away the magnetic flux so that the QED equation, divB = 0, is 
satisfied identically. The concrete form of the 'string' depends on the concrete way 
one parametrizes the continuation of the physical vector u^{ti,t 2 ) to the unphysical 
dimension labelled by a. Since the action is given by the 'winding number' eq. ([22]), 
different parametrizations may differ only by multiples of 2ir. Let us stress that the 
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correct coefficient in ([27]) representing the Dirac quantization condition, follows from 
the normalization factor 1/247T 2 of the winding number, see eq. (22). 



4 Non-abelian monopoles 

Writing down the charge-monopole interaction in the Wess-Zumino form we have ac- 
tually integrated off the gauge fields, leaving only the charge and monopole trajectories 
as the dynamical variables of the theory. Therefore, exp(£<Sj n t) is, in fact, what is called 
the generating functional of the theory, depending on the external charge current j^(x) 
determined by the charge trajectory x^(r). In case of the Yang-Mills theory the exter- 
nal colour current j^(x) is an adjoint matrix. Some time ago we have shown Q that 
the Yang-Mills generating functional has specific properties which follow from gauge 
invariance. It can be written as a sum of two pieces: one is gauge invariant (call it Wi) 
and can depend only on the diagonal part of the current d^, 

W = S(T)%(T)£ft(T) (28) 

(one can always locally rotate the adjoint matrix to make it diagonal), while the second 
piece is gauge non- invariant, and depends on the 'angle' variables S(t). According to 
ref. 11 the general form of the generating functional in the Yang-Mills theory for a 
closed loop of a point-like current can be written as 

W%}=W 1 [dl(r)] + J J drSS\ (29) 

In calculating the vacuum average of the Wilson loop one actually has to substitute 
the appropriate colour current j^, as induced by the specific loop considered, into the 
generating functional of the theory. The appropiate current of the Wilson loop can be 
read off from (|]) (or eq. (J3j) for a general gauge group). It follows from eq. (Q) that 
the diagonal part of the probe current is nothing but 

d*(x) = 2J Jdr x^r) 6^(x a - x a (r)), (30) 

where x a {r) is the trajectory of the colour charge, while its 'angle' part S{t) is nothing 
but the gauge transformation matrix over which one has to integrate in eq. (Q). Com- 
paring eqs.(|4,29) one observes that the second (gauge non-invariant) pieces are exactly 



cancelled. This has been the conclusion of ref. Q. Therefore, the integration over 
all possible gauge transformations S(t) in eq. (0) becomes trivial, so that one is left 
only with the gauge-invariant piece of the generating functional, W\. If one assumes 
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that the dynamics of large Wilson loops is governed by monopoles, this W\ is just the 
Wess-Zumino-type charge-monopole interaction term, eq. (p2|), where one has to use 
the diagonal part of the colour charge current, given by eq. (|30|). 

It is less evident how to introduce colour monopoles, however one can think of 
proceeding in a similar way as with the colour charges. Namely, the interaction of 
colour magnetic charges with the dual Yang-Mills potential is given by a dual 
Wilson loop (cf. eq. (||)), 



L m = J DS{r)exp iK J dr Tr r 3 (SBS j + iS&) 



(31) 



where B(t) is the tangent component of the dual Yang-Mills potential along the 
monopole loop, S(t) is the unitary matrix of the instant colour orientation of the 
monopole, K = ^, 1 is the colour 'spin' of the monopole. As in the case of the 
similar eq. (Q) one has to integrate actually over unit orientation vectors m(r) from 
the coset space, rather than over the group elements S(t). 

One can introduce the field theory of monopoles in the path-integral formalism 
and write the partition function as a sum over arbitrary numbers of monopole loop 
trajectories, with their colour orientation vectors m(r) living on those trajectories: 

Z K = E H ft / Dm ^ T ) I D y^ T ) MvUr)l (32) 

n n - i J J 

where A^-fy^r)] are some weight functions. For the free particles of mass mx the 
weight function is Ak(1) = exp(— m^l) where I = / dr^yf a is the monopole loop 
length. With such a weight functional eq. ([32]) is equivalent to the partition function 
of the field theory of free scalar Yang-Mills particles with mass txik and colour 'spin' 
belonging to a certain representation labelled by K. One can take into account that 
monopoles can, in principle, belong to different representations of the gauge group, 
however for simplicity we consider only some particular K. In a theory of interacting 
monopoles the weight functional is, of course, more complicated. 

The main idea of this paper is to avoid the notion of monopole fields and of dual 
Yang-Mills potentials (as well as the usual ones), and to work directly with charge 
and monopole trajectories. The Wess-Zumino-type interaction term (^) solves the 
problem of eliminating the dual as well as the usual Yang-Mills potential, therefore 
it can be understood as the generating functional - not only for the electric charges 
but also for the magnetic ones, since they enter the interaction term on equal footing. 
Therefore, everything said above about the non-abelian electric currents applies to the 
non-abelian magnetic currents. The Wess-Zumino term containing the instant colour 
orientation of monopoles (the second term in eq. (|3l|)) has to be cancelled exactly by 
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the gauge non- invariant term of the generating functional similar to that of eq. (p9|). 
Then the integration over rrij in the monopole partition function (^) becomes trivial, 
and one is left with the abelian charge-monopole interaction term where the abelian 



currents are given by eq. (30) for the electric charges and by 



<%(x) = 2K J dr y^r) fiW(x a - y a (r)), (33) 

for the magnetic currents, where y a (r) are the trajectories of monopoles. 

We arrive thus to a simple recipe to compute the non-abelian electric Wilson loop 
in the background of fluctuating quantum monopole fields - without introducing such 
fields at all. It is given by the exponent of the charge-monopole interaction (p^), 
averaged with the monopole partition function (|32|). We obtain for the averaged 
Wilson loop whose trajectory is denoted by Xu,(t): 

1 1 n r 

( L Mt))) = ^E^Il / %(^)) A K [y ia (r)] 

n ' i 

x e^ V {{i(2J){2K)S int {x^y ilx }} . (34) 

The generalization to arbitrary gauge groups and representations of both charges 
and monopoles is straightforward, according to eq. (||): one has to replace 2Jr 3 (and 
2Kt 5 ) by J2 m iHi where Hi are r generators from the Cartan subalgebra and m is the 
r-dimensional eldest weight of the representation to which the electric charge belongs 
(and similarly for the magnetic charge). 

Eq. (S3) relates directly the Wilson loop to the integral over monopole paths. If 
confinement, i.e. the area behaviour of the Wilson loop, is due to the monopole con- 
densation, it should be seen from this formula. Both monopole condensation and 
confinement can be thus formulated in completely gauge-invariant terms: it is a state- 
ment about the weight functional AK\yi a {j)\. For example, one can easily check that 
the weight Ak ~ exp(— ml) corresponding to monopoles being free massive particles, 
does not lead to the area law: monopoles should be at least effectively massless. Long 



loops are absolutely necessary to derive the area law from eq. (34) [14]. 



5 Conclusions 

We have derived a new non-abelian Stokes theorem for the Wilson loop. The path 
ordering of the Wilson loop is exchanged for the integration over an auxiliary scalar 
field from the coset SU(N)/[U(1)]^ N ~ 1 ^ space; its meaning is the instant direction of 
the colour charge. With the help of that field one can construct a gauge-invariant field 
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strength to which the Stokes theorem applies. The construction needs an introduction 
of a non-uniquely defined Wess-Zumino term for the auxiliary field, which fixes the 
group representation of the Wilson loop considered. 

We have presented the interaction between quark and monopole in a Lorentz- 
invariant, gauge-invariant and reparametrization-invariant form. It is also given by 
a Wess-Zumino-type action; its non-uniqueness is due to the arbitrariness of how one 
draws the Dirac string. However, under the change of the string direction it can only 
produce phase factors exp(27rm) = 1. 

The average of the quark Wilson loop over the vacuum filled with fluctuating 
monopoles has been related directly to path integrals over monopole trajectories. This 
enables one to formulate in a gauge invariant way - in the language of monopole world 
lines - what precisely does the alleged monopole condensation mean, and how does 
one get the area behaviour of the Wilson loop out of that [14]. 
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